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We reply to the Comment by Achim Rosch 1) who challenges our finding in Ref. 2) with 
regard to the validity of the Luttinger theorem in the cases of Mott insulating A-particle 
ground states (even) when the chemical potential used in applying this theorem coincides 
with the zero-temperature limit of the chemical potential satisfying the equation of state 
corresponding to N particles. Rosch further argues that the strong-coupling expression for 
the single-particle Green function presented in Ref. 3) and analyzed in Ref. 2) does not 
imply destruction of the Mott insulating state at half-filling as a result of an arbitrary weak 
hopping contribution that breaks particle-hole symmetry and therefore suggests that our 
conclusion in Ref. 2) to the contrary were incorrect. Here we show the shortcomings of 
Rosch's arguments. 
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In Ref. 2) (Sec. 6.1 herein) we explicitly considered two specific cases which 
had been investigated by Rosch in Ref. 3): (1) the local or zero-hopping limit, and 
(2) the non-local case corresponding to small hopping energy dispersion r^, Vfc, in 
comparison with the effective interaction energy U. The reader is referred to Ref. 3) 
for the details of the model to which these cases correspond. 

In the local limit, for the Green function*) G(k;z) and self-energy S{k\z) one 
has 2 )' 3 ) 

£ tac(2) 4+<^. (2) 

Observing that 2 )' 3 ) (below, the last equality applies only for < U /2): 

dz ~ 8 ~ m+ioo ^ z _ q _ i 

^— Gi oc (z) — E loc (z) = - — Ui oc (z)— Gioc(z) = -sgn(/x), (3) 

and given that at zero temperature for fj, inside the Mott gap (—U/2,U/2) the 
number of particles is independent of the value of fi, Rosch 3 ) correctly 2 ) concluded 
that the Luttinger- Ward identity, 4 ) and thus the Luttinger theorem, 4 )**) breaks down 
in the local limit. Although Rosch conceded that the Luttinger theorem applies 



*- ) Throughout we identify h and k B with unity. As a result of fc B = 1, in the following j3 = 1/T. 
**' In Ref. 2) (Sec. 4.3 herein) we demonstrated that the Luttinger theorem applies if and only if 
the Luttinger- Ward identity is valid, and that deviation from this result implies that the underlying 
ground state is pathological. 



2 



Behnam Farid 



for jjL = fj, h (in the local limit, /x L = 0), nonetheless this does not distract from 
the fact that for insulating ground states no specific aspect in the proof of the 
Luttinger theorem seems to depend on one distinct value of \x so that the validity 
of the Luttinger theorem in the case at hand for // = and its failure for all other 
fi £ {—U /2, U /2) may be indeed characterised as failure of this theorem. 2 ) 
With 

( m = H + iuj m , (4) 

where io m , m = 0, ±1, . . . , is the mth Matsubara frequency, in Ref. 2) we obtained 
that 

\ E GiocW = \ [2 tanh (St) - tanh - |)) - tanh + |)) 

= $(J3H,/3U). (5) 

The relationship between the sum on the left-hand side (LHS) of Eq. (5) and the 
integral on the LHS of Eq. (3) is established through the equality 



-i rfi+ioc J /- 

\im-J2MCm)= ^lim/zKC). ( 6 ) 



This correspondence, which plays a significant role in the derivation of the Luttinger- 
Ward identity, 4 ^ follows from the consideration that d£ m = ( m +i — Cm = 2m/ j3 — > 
for (3 — ► oo, so that in this limit the sum on the LHS of Eq. (6) is the Riemann sum 
of the integral on the right-hand side (RHS). 

Although the Green function and the self-energy on the LHS of Eq. (5) are the 
zero-temperature limits of their corresponding finite-temperature counterparts,*) we 
showed in Ref. 2) that as (3 — > oo, for /x inside (-U /2, U /2) and sufficiently far from 
±U/2, to exponential accuracy the RHS of Eq. (5) coincides with the expression that 
one would obtain on employing in the LHS of this equation the finite-temperature 
counterparts of G\ oc (z) and U\ oc (z). In fact, for (3 — ► oo to the same accuracy one 
has 2 ) 

\ E gioc(Cm) m :f m) ~ <fr(p» = \ tanh A) for - oo. (7) 

The combination /3 x /i on the RHS of this expression makes explicit that the zero- 
temperature limit of the sum on the LHS is approached when S> 1/fx. Conse- 
quently, for (i — ► the low-temperature regime corresponds to increasingly larger 
values of (3. Aside from this fact, one notes that since 

lim lim = --. lim lim &(/3fj,) = -, lim lim <F(/3//) = 0, (8) 

the validity or failure of the Luttinger-Ward identity, and thus of the Luttinger theo- 
rem, crucially depends on the way in which the zero-temperature limit is effected. In 



*' With reference to Eq. (6), the LHS of Eq. (5) corresponds to lim^oo \ ~^2 m faa{C,n 
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Ref. 2) we referred to the first two repeated limits (§§ 302-306 in Ref. 5)) in Eq. (8) as 
'false', or 'spurious', limits.*) The last repeated limit in Eq. (8) in conjunction with 
the expression in Eq. (7) establish that in the event that the zero-temperature limit 
of the chemical potential satisfying the equation of state at finite temperatures (i.e. 
fi/3 = /i(/3, N, V), where V is the volume of the systems in the grand-canonical ensem- 
ble under investigation) is vanishing, the Luttinger- Ward identity, and therefore the 
Luttinger theorem, applies in the local limit on identifying fx with fioo = lim^oo (ip 
prior to effecting the zero-temperature limit. 

In Ref. 2) (Sec. 6.1.3 herein) we explicitly demonstrated that in the local limit 
fj,p approaches zero faster than 1/(5 for (3 — > oo. In consequence of this, not only 

lim lim #(/fyt) = lim ^(P^oo) = 0, (9) 

but also 

lim &(Pw3) = 0. (10) 

f3—*oo 

We have thus explicitly demonstrated 2 ) that, in the local limit the Luttinger- Ward 
identity, and thus the Luttinger theorem, is obtained on equating /i with either 
/ioo or fj,p prior to evaluating the zero-temperature limit. This result unequivocally 
demonstrates that in the local limit the genuine breakdown of the Luttinger theorem 
is averted by identifying the value of the thermodynamic variable \x with //oo- 

In 

Ref. 2) we present arguments showing that the validity of the Luttinger theorem for 
A* = A^oo is general and not restricted to the local limit considered here. 

As for the non-local cases, to linear order in t^/U for the Green function corre- 
sponding to these, Rosch 3 ) obtained that 

G(k; z) = = = . (11) 

z + U/2 - T fc - S\ oc {z) 

With 

oo±(k): = — ±-(t% + U 2 ) 1 ^ 2 = — ± — + oCi) as (12) 

denoting the poles of the Green function in Eq. (11), for the single-particle spectral 
function A(k;u>) corresponding to this Green function one has 2 ) 

A{k;u) = LzWM±gg! , ( „ _ a _ m + ^Mig , ( „ _ w+m 

(13) 

It is to be noted that the sum of the weights of the two delta functions in this 
expression is equal to unity, as befits a properly normalised single-particle spectral 
function. 

For a given value of the chemical potential fx, at zero temperature for the mean 
value of the number of particles per site per spin spices, n, one has 

n = wT, T ^A(k;co), (14) 

k 00 



*' Here, as in Ref. 2), the terms 'false' and 'spurious' are used in their technical sense and do 
not connote false or spurious mathematical operations carried out by a particular researcher. 
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where Af h is the number of lattice sites and the factor 2 accounts for the two orbitals 
per site in the model under consideration. Let now 

J?_:= — minife > 0, i? + :=maxTfc > 0. (15) 
fc fc 

Introducing the normalised density-of-states function 

fc 

from the expressions in Eqs. (13) and (14) one trivially obtains that 2 ) 
n = 1 - C for maxw-(fe) < n < mincj+(fc) <==> ~ + ^± < fi < ^ - (17) 



where 



We note that 2 ) 



fc v y r fc^ v/ 2 2 ~ r ~ 2 2 



do, ^ , . (18) 

✓ O . T-Tl\1 /O \ ' 



+ [/2)l/2 



/" dw = 1, y dw w = 0. (19) 

For symmetric cases, corresponding to T>{— u) = T>(u), Vw, the constant C is 
identically vanishing. For asymmetric cases and U > max(f?_, ft + ) one has the 
following uniformly convergent series 2 ^ 



C = Xi \- X 3 + -4- X 5 - . . . , where J ? := / dcu V(u) J . (20) 

U 2U 3 8U 5 



On account of the second expression in Eq. (19) one has X\ = 0, so that C is at the 
largest of the order of (loq/U) 3 where ojq G (—!?_, ,17 + ); for convenience, below we 
shall refer to C as being at the largest of the order of (rfc/C/) 3 . Note that o C / 
implies a net amount of total spectral weight transfer from the band uj^(k) to the 
band u>±(k) when C ^ 0. 

Since for asymmetric cases C is not identically vanishing (it may however be 
vanishing for some 2 ) discrete values of U), from Eq. (17) it follows that for (i in the 
interval given in Eq. (17) the ground state of the system under consideration cannot 
be half-filled, i.e. n = 1 is ruled out; the requirement for the ground state to be 
half-filled, turns the ground state metallic. 2 ) 

Rosch 1 ) considers the conclusions in the previous paragraph as incorrect. He 
argues that the Green function presented in Eq. (11) being correct only to linear 
order in Tk/U, C, which at the largest is of the order of (rk/U) 3 (see above), should 
be identified with zero. We, on the other hand, reason that the half-filled metallic 
ground states predicted by the Green function in Eq. (11) being induced by merely 
C / 0, irrespective of how small \C\ may be, signifies the fact that these half-filled 
metallic ground states are not to be taken as physically viable. In the closing part 
of Sec. 6.1.5 in Ref. 2) we stated that: 
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"Summarising, since not until the third order in Tk/U can C deviate from zero in 
asymmetric cases, we have thus the clearest evidence that first-order results are in 
principle incapable of establishing break-down of the Luttinger theorem in asymmetric 
cases." 

In Ref. 2) (the closing paragraph of Sec. 6.1.13 herein) we further acknowledged the 
observation by Rosch 3 ) that for [i in the vicinity of =Lt//2 the leading-order pertur- 
bation expansion resulting in the Green function in Eq. (11) becomes inadequate. 

Now let us investigate the consequences of Rosch's proposal of identifying C with 
zero. Since (cf. Eq. (20)) 

r k l/^3 + 3 ( 3 )5 _ for W <1; (21) 



{t 2 + U 2 ) 1 / 2 U 2 y u' % K U J U 

it follows that identification of C with zero amounts to approximating Tfc/(r| + 
U 2 ) 1 / 2 by the leading-order term Tk/U. Following Eq. (13), this approximation 
combined with the concomitant leading-order approximations of u>±(k) result in the 
approximate single-particle spectral function 

A'(k; u) = 1 ~ T 2 k/U S(u - r k /2 + U/2) + 1 + r ^ U S (u - r k /2 - U/2). (22) 

One readily verifies that this function corresponds to the approximate single-particle 
Green function 

G'(k; z) = -( ^- + = 1 , . (23) 

2\ z -T k /2 + U/2 z- Tk /2-U/2J z - r k + [r 2 - U 2 )/(Az) 

With Go(k; z) = l/(z — T k + U /2), from the Dyson equation one immediately infers 
that G'(k;z) corresponds to the self-energy 

E'(k;z) = E loc (z)-^. (24) 

In contrast to G(k; z), G'(k; z) is thus seen to correspond to a non-local self-energy. 
It may be instructive to express r\j z as {r k jz)T k whereby one may view r\j z as 
'linear' in r k with a diverging coefficient for z — > 0. This observation is noteworthy 
in that the leading-order hopping contribution to G\ oc (z), resulting in the Green 
function in Eq. (11), is the — T k encountered in the denominator of the expression 
on the RHS of this equation. 

Not surprisingly, one can explicitly show that the half-filled ground state corre- 
sponding to G'(k; z) is Mott insulating, exactly as is the case for G\ oc (z) and IJ\ oc (z), 
and indeed for Gk; z) and X!\ oc (z) when T>(uj) is symmetric. 2 ) This interesting obser- 
vation notwithstanding, it should be evident that the considerations in Ref. 2) have 
no bearing on G'(k; z), but on the G(k; z) presented in Eq. (11). 

We have thus unequivocally established that a theoretical framework in which 
C is identified with zero, on account of it being of higher order than linear in the 
small parameter Tk/U, is strictly distinct from the framework in which the G(k;z) 
in Eq. (11) is the exact Green function. In this light, we believe that Rosch 1 ) misin- 
terprets a number of matters by stating that 1 - 1 



6 



Behnam Farid 



"His [Farid's] argument is based on a surprising result of his calculations: he claims 
[3] that an arbitrary small breaking of particle-hole symmetry transforms the half- 
filled Mott insulator into a metal, or, equivalcntly, that the particle-hole asymmetric 
system in not half-filled if the chemical potential is located within the gap! In our 
opinion, this is obviously wrong. For example, it contradicts the observation that 
small perturbations have no effects in systems with a finite gap (in the two-band Mott 
insulator under consideration both the charge and the spin gap are finite). In the 
appendix we sketch the formal argument which can be used to prove this." 
The "surprising" result to which Rosch 1 ) refers is presented in Eq. (17) above, which 
arises from a rigorous treatment of the Green function in Eq. (11). To contrast 
this result, corresponding to an approximate Green function, with one deduced from 
exact considerations, and subsequently characterise the calculations leading to the 
"surprising" result as incorrect is not only logically unsound, but is grossly unfair. 
Nowhere in Ref. 2) have we indicated or implied that the conclusions that Rosch 
identifies as "surprising" and "not correct" should be viewed as exact in the absolute 
sense; the validity of these conclusions is relative to a framework in which the Green 
function in Eq. (11) is held as exact. 

Similarly, the statement with regard to the stability of the half-filled Mott insu- 
lating phase against small asymmetric perturbations cannot have any place within 
the framework where the Green function in Eq. (11) is exact. The exact considera- 
tions presented in the appendix to Ref. 1), which are similar to those explicating the 
robustness of the Mott phase against small-amplitude hopping terms in the Bose- 
Hubbard Hamiltonian (§ 10.1 in Ref. 6)), does not apply here; such considerations 
in principle apply only if the Green function under consideration is the exact Green 
function. Stated differently, with H denoting the Hamiltonian of the system under 
investigation, it is required that G(k; z) correspond to the exact iV-particle ground 
state of H. As it stands, the G{k;z) in Eq. (11) not only does not correspond to 
the A^-particle ground state of H, with absolute certainty it does not correspond 
to the iV-particle ground state of any interacting Hamiltonian operating in infinite- 
dimensional N- and N ± 1-particle Hilbert spaces. 

To summarise, those conclusions arrived at in Ref. 2) with which Rosch 1 ) specif- 
ically disagrees, are exact. The source of Rosch's contention appears to be two-fold. 
Firstly, Rosch believes, unwarrantedly as our above calculations show, that the con- 
stant C should be identified with zero. Secondly, for reasons that are not apparent 
to us, Rosch 1 ) asserts that the "surprising" result in Eq. (17) is not to be trusted 
on account of it contradicting the robustness of the Mott insulating state against 
small hopping terms in the underlying Hamiltonian. As for the first point raised by 
Rosch, we have explicitly shown that identification of C with zero leads to trans- 
formation of the G(k; z) in Eq. (11) into an entirely different Green function, i.e. 
the G'(k;z) in Eq. (23); evidently, the considerations in Ref. 2) have no bearing, 
whatever, on the properties of G'(k;z). As for Rosch's second point, we have em- 
phasised that our "surprising" finding is only so by assuming that the Green function 
in Eq. (11) were exact in the absolute sense. We have further underlined the fact 
that the considerations in the appendix to Ref. 1) would apply if the Green function 
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in Eq. (11) corresponded to the exact ground state of the interacting Hamiltonian 
under considerations. 

I should like to thank A. Rosch for kindly providing me with a draft of his 
Comment in advance of its publication. 
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